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Abstract 



In this paper, a new discrete scheme for Gaussian curvature is pre- 
^D \ sented. We show that this new scheme converges at the regular vertex 

with valence not less than 5. By constructing a counterexample, we 
also show that it is impossible for building a discrete scheme for Gaus- 



sian curvature which converges over the regular vertex with valence 



<^ ' 4. Moreover, the convergence property of a modified discrete scheme 

for the Gaussian curvature on certain meshes is presented. Finally, 
asymptotic errors of several discrete schemes for Gaussian curvature 



are compared. 

AMS Subject Classifications: Primary 68U07, 68U05, 65S05, 53A40. 
Keywords: Discrete Gaussian curvature, discrete mean curvature, geometric 
modeling. 



1 Introduction 

Some applications from computer vision, computer graphics, geometric mod- 
eling and computer aided design require estimating intrinsic geometric invari- 
ants. It is well known that Gaussian curvature is one of the most essential 
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geometric invariants for surfaces. However, in the classical differential ge- 
ometry, this invariant is well defined only for C 2 smooth surfaces. In mod- 
ern computer-related geometry fields, one often uses C° continuous discrete 
triangular meshes to represent smooth surfaces approximately. Hence, the 
problem of estimating accurately Gaussian curvature for triangular meshes 
is raised naturally. 

In the past years, a wealth of different estimations have been proposed 
in the vast literature of applied geometry. These methods for estimating 
Gaussian curvature can be divided into two classes. The first class is based 
on the local fitting or interpolation technique [2j [5j [6j [7J [20] , while the second 
class is based on discretization formulations which represent the information 
about the Gaussian curvature [H El El 021 EH EB] • In this paper, our focus is 
on the methods in the second class. The main aim of the paper is to present 
a new discrete scheme for the discrete Gaussian curvature which converges 
at the regular vertex with valence not less than 5. 

Let M be a triangulation of smooth surface S in M. 3 . For a vertex p of 
M, suppose {pi}" =1 is the set of the one-ring neighbor vertices of p. The set 
{PiPPi+i} (i = 1,- ■ ■ ,n) of n Euclidean triangles forms a piecewise linear 
approximation of S around p. Throughout the paper, we use the following 
conventions p n +i = Pi and p = p n - Let ji denote the angle Zpjpp i+1 and 
the angular defect at p be 2n — £\ 7^. 

A popular discrete scheme for computing Gaussian curvature is in the 
form of # 7 \ where E is a geometry quantity. In general, one selects E 
as A(p)/3 and obtain the following approximation 

C (i) ._ 3(27r-£ t 7*) 

A(p) ' [L) 

where A(jp) is the sum of the areas of triangles pjppj + i. In [Jj, another scheme 

G (2) ;= 2lT-£ <7i (2) 



Op 



is given, where 



s ? ■= E z^r [ww " ^P^ 2 + vl 
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4-1; 



— ' 4 sin 7j 
is called the module of the mesh at p. In [16], the discrete approximation 



G^ is modified as 

G 0) := 27r-^ i7i ^ 

~ £). area(p 4 pp i+ i) - | £\ cot(7*)df ' 

where di is the length of edges PiPi+i. There are several different points of 
view for explaining the reason why the angular defect is closely related to 
the Gaussian curvature with including the viewpoints of Gaussian-Bonnet 
theorem, Gaussian map and Legendre's formula (see the next section for 
details). 

Asymptotic analyses for the discrete schemes have been given in [TJ [7J 
ITS] . In [7], the authors show that for the non-uniform data, the discrete 
scheme G^ does not convergent to true Gaussian curvature always. In [JJ, 
Borrelli et al. prove that the angular defect is asymptotically equivalent to a 
homogeneous polynomial of degree two in the principal curvatures and show 
that if p is a regular vertex with valence six, then the scheme G^ converges 
to the exact Gaussian curvature in a linear rate. Moreover, Borrelli et al. 
show that 4 is the only value of the valence such that the angular defect 
depends upon the principal directions. In [TS], Xu proves that the discrete 
scheme G^ has quadratic convergence rate if the mesh satisfies the so-called 
parallelogram criterion, which requires valence 6. Therefore, one hopes to 
construct a discrete scheme which converges over any discrete mesh. But in 
|21j . Xu et al. show that it is impossible to construct a discrete scheme which 
is convergent for any discrete mesh. Hence, we have to be content with the 
discrete schemes which converge under some conditions. According to the 
past experience [H [SJ [2TJ , we regard a discrete scheme desirable if it has the 
following properties 

1. It converges at a regular vertex, at least for sufficiently large valence 
(the definition of the regular vertex will be given in Section 2); 

2. It converges at the umbilical point, i.e., the points satisfying k m = ku 
where k m and k M are two principal curvatures. 

As stated before, the previous discrete schemes, including G^ l \G^ 2 ' and 
G^\ only converge at the regular vertex with valence 6. In [1], a method for 
computing the Gaussian curvature at the regular vertex with valence unequal 
to 4 is described. But the method requires two meshes with valences n\ and 
n 2 { n i 7^ 4, n 2 7^ 4, rii 7^ n 2)- I n this paper, we will construct a discrete 
scheme which converges at the regular vertex with valence not less than 5, 



and also at the umbilical points with any valence. Moreover, the discrete 
scheme requires only a single mesh. Hence, the new discrete scheme is more 
desirable. Furthermore, we show that it is impossible to construct a discrete 
scheme which is convergent at the regular vertex with valence 4. Therefore, 
for the regular vertex, the convergence problem remains open for the case 
of vertex with valence 3. We also want to point out that we discuss the 
pointwise convergence in this paper, which is different with the convergence 
in norm as discussed in [H [10] . 

The rest of the paper is organized as follows. Section 2 describes some 
notations and definitions and Section 3 shows three different viewpoints for 
expressing the relation between the angular defect and Gaussian curvature. 
In Section 4, we study the convergence property of the modified discrete 
Gaussian curvature scheme. We present in Section 5 a new discrete scheme 
and prove that the scheme has good convergence property. In Section 6, for 
the regular vertex with valence 4, we show that it is impossible to build a 
discrete scheme which is convergent to the real Gaussian curvature. Some 
numerical results are given in Section 7. 

2 Preliminaries 




Fig.l. Notations. 



In this section, we introduce some notations and definitions used through- 
out the paper (see Fig. 1). Let S be a given smooth surface and p be a point 



over S. Suppose the set {pjPPj + i}, i = 1, • • • ,n, of n Euclidean triangles 
form a piecewise linear approximation of S around p. The vector from p to 
Pi is denoted as ppj. The normal vector and tangent plane of S at the point 
p is denoted by n and II, respectively. We denote the projection of p» onto 
II as qj, and define the plane containing n, p and p^ as IT. Then we let Ki 
denote the curvature of the plane curve S D H at p. The distances from p 
to pj and q« are denoted as rji and l iy respectively. The angles Zpjpp i+1 and 
ZqiPQi+i are denoted as 7; and /?$, respectively. The two principal curva- 
tures at p are denoted as k m and Um- Let r\ = maxj7?j. The following results 
are presented in pfl [8j [18] : 

^- = 1 + 0(77), fr = % + 0( V 2 ), (4) 

Vi 



E 



WiPPi 



J2^ + o( v % (5) 



where Wi 6 R. 

Now we give the definition of the regular vertex using the notations in- 
troduced above. 

Definition 2.1. Let p be a point of a smooth surface S and let Pi,i = 
1, • • • , n be its one ring neighbors. The point p is called a regular vertex if 
it satisfies the following conditions 

(1) the A = f , 

(2) the rjiS all take the same value 77. 

Remark 2.2. We can replace (1) in Definition 12.11 by requiring the 7$ all 
take the same value. Since A = 7« + 0(r] 2 ), all the results in the paper hold 
also for the alternative definition. 

3 Angular Defect and Gaussian Curvature 

In this section, we summarize three different viewpoints for expressing the 
relation between angular defect and Gaussian curvature. These viewpoints 
have been described in different literature [TJ, [16l HE]- We collect them there. 
Throughout the section, we use G^ (p) to denote the discrete Gaussian cur- 
vature at p, which is obtained using G^ l \ 



3.1 Gaussian- Bonnet theorem viewpoint 

Let Dbea region of surface S, whose boundary consists of piecewise smooth 
curves I^s. Then the local Gaussian-Bonnet theorem is as follows 

// G{p)dA + Y] / k g {T j )ds + y2a j = 2'ir, 

JJd j Jv 3 . 

where G{j>) is the Gaussian curvature at p, k g (Tj) is the geodesic curvature 
of the boundary curve 1^- and atj is the exterior angle at the jth corner point 
Pj of the boundary. If all the TjS are the geodesic curves, the above formula 
reduces to 

jf G{p)dA = 2vr - J2 a r ( 6 ) 

Let M be a triangulation of surface S. For vertex p of valence n, each triangle 
P*PPi+i can be partitioned into three equal parts, one corresponding to each 
of its vertices. We let D be the union of the part corresponding to p of 
triangles PiPP i+ i- Note that Y^ili = J2j a j- Assuming G(p) is a constant on 
D, and using flS]), we have G(p) can be approximated by G^(p). 

3.2 Spherical image viewpoint 

We now introduce another definition of Gaussian curvature. Let D be a 
small patch of area A including point p on the surface S. There will be a 
corresponding patch of area / on the Gaussian map. Gaussian curvature at 
p is the limit of ratio lim^^o 4- 

Let us consider a discrete version of the definition. The Gaussian map im- 
age, i.e. the spherical image, of the triangle PiPP i+ i is the point |i( P I P ')x( - P ')|i 
Join these points by great circle forming a spherical polygon on the unit 
sphere. The area of this spherical polygon is lix — J^ 7$. Same as the above, 
each triangle is partitioned into three parts, one corresponding to each vertex. 
Then the Gaussian curvature can be approximated by (?W(p). 

3.3 Geodesic triangles viewpoint 

Let T = ABC be a geodesic triangle on the surface S with angles a,/3, 7 
and geodesic edge lengths a, b, c. Let A'B'C be a corresponding Euclidean 



triangle with edge lengths a, b, c and angles a', j3', 7'. Legendre presents the 
following formulation 

a - a' = area(T)— \-^- + o(a 2 + b 2 + c 2 ), 

where area(T) is the area of the geodesic triangle ABC, G(A) is the Gaussian 
curvature at A. 

Using Legendre's formulation for each triangles with p as a vertex, we 
arrive at the estimating formula G^{p) again. 

4 Convergence of Angular Defect Schemes 

In [18], Xu gives an analysis about the scheme G^ and proves that the 
scheme converges at the vertex satisfying so-called parallelogram criterion. 
A numerical test shows that the scheme does not converge at the regular 
vertex with valence unequal to 6 and at umbilical points. In [I], Borrelli et. 
al. give an elegant analysis about the angular defect. They show that if the 
vertex p is regular, then the angular deficit is asymptotically equivalent to 
a homogeneous polynomial of degree two in the principal curvatures with 
closed form coefficients. Moreover, they present another angular scheme 
q{2) ._ j^ 7 ' . In fact, using the law of cosine, we have 

Dp 



- Y^ area( Pi pp 4+1 ) - - ^ cot( 7i )< 2 



E 



1 . 1 COS 7j . o 2 . 

-ViVi+i sm7i - 77. + T] i+1 - 2r]iT] i+1 cos 7^ 

4 8 sin 7j 



El r cos7i 2 2 ■ 
-T-. ViVi+i TTKHi + Vi+i) 
4sm7j L 2 . 



b p . 



This shows that G^ and G^ are equivalent, which means these two schemes 
obtain the same value for the same triangular mesh. 

In [T8], the author proves that the discrete scheme G^ l > has quadratic 
convergence rate under the parallelogram criterion. In the following theorem, 
we shall show that the discrete scheme G^ has also quadratic convergence 
rate under the same criterion. 



Theorem 4.1. Let p be a vertex of M with valence six, and let Pj,j = 
1, • • • , 6 be its neighbor vertices. Suppose p and Pj,j = 1, • • • ,6 are on a 
sufficiently smooth parametric surface F(£i,£2) £ K 3 , and there exist u,Uj 6 
M 2 such that 

p = F(u), pj = F(Uj) and U, — u = (uj_i — u) + (iij+i — u), j — 1, • • • , 6. 

T/ien 

"" ' ; G( P ) + 0(r 2 ), 



|A(p,r)-jE,oot( 7 i(r))d?(r) 

where, G(p) is the real Gaussian curvature o/F(u) at p, 

^(P, r ) : = 5Z area [Pi( r )PPi+i( r )]' Pi( r ) := F K( r ))> 

i 

and Uj(r) = u + r(uj — u), z = 1, • • • ,6. 

Proof. Let 

A(p, r) = a r 2 + a x r 3 + 0(r 4 ) (7) 

and 

^f^ " I E cot( 7 ,(r))^ 2 (r) = 6 r 2 + ^r 3 + 0(r 4 ) 

i 

be the Taylor expansions with respect to r. According to Theorem 4.1 in 

3^ = G(p) + 0(r»). 

Hence, to prove the theorem, we need to show b = a /3, b\ = ai/3. Ac- 
cording to [12], we have Oi = 0, which implies that we only need to prove 
b = ao/3,6i = 0. 

Note that the u, u.j,j = 1, ... 6, satisfy the parallelogram criterion. With- 
out loss of generality, we may assume u = [0, 0] T , ui = [1, 0] T . Then there 
exists a constant a > and an angle 9 such that 

u 2 = [a cos 9, a sin 9] T . 

Hence, u 3 = [a cos 9 — 1, a sin 9} T , Uj +3 = — Uj, j = 1, 2, 3. Let 

U J = s J d J = s A9j, Ijf, 3 = 1, • • • , 6, 



where Sj = lluJI and ||cL|| = 1. Then, we have 



s 1 — 1, s 2 — a, s 3 = Va 2 - 2ac+ 1, s 4 = s 1 , s 5 = s 2 , s 6 = s 3 , 

9\ = 1, #2 = c, 5-3 = (oc - l)/s 3 , 5-4 = ~9l, 95 = ~92, 96 = ~93, 

h — 0, I2 — t, I3 = at/ S3, U = — /1, I5 = — 12, l& = —I3, 
where (c, t) := (cos 9, sin 0). Note that 

1 6 / 

A (P> r ) = 2EV IIPi( r )-Pll 2 HPi+i( r )-Pll 2 - (Pi( r )-P>Pj+i( r )-p) 2 > 

cot( 7 ,(r)) = ( Pj (r)-p,p J+1 (r)-p) 

VllPi( r ) -p|| 2 IIPj+i( r )-p|| 2 - (Pj( r ) -P>Pi+i( r ) -P) 2 

(9) 
d%r) = || Pj (r) - p|| 2 + ||p j+ i(r) - p|| 2 - 2( Pj (r) - p, p i+1 (r) - p). (10) 



Let F j denote the kth order directional derivative of F in the direction 
a j 

dj. Then using Taylor expansion with respect to r, we have 

l|Pi(0 - P.ll 2 = sy(F dj ,F dj ) + ^(F^F 2 ,) + ± 8 y {¥%,¥%) 



+- s y(F dj ,F%) + - s y(Fi,F%) + - s y(F dj ,F%) + oy), (n) 



± s y { F dj ,Fi) + lsy(Fi,Fi) + ± 2 

and 



(Pj(r)-p,p i+ i(r)-p) 

s r s J+1 r 2 (F dj ,F d3+1 ) + -s 3 s) + y{F dj ,F\ +i ) + -s 2 jSj+1 r 3 {F dj+i: F 2 .) 

k^V^^F 2 ,) + ^V 4 (F d ,,F|. +i > + i^ 3 ^ + ir 4 (F, j+1 ,F d3+1 ) 






+ ^U a .ir S (F i ,,F%J + ±s J+1 sy(F> t ,F Ji J + 0(r<>). (12) 

To compute all the inner products in the two equations above, we let 
_ aF(ft,6,) _ g 2 F(ei,6) . d 3 F d 4 F 

% ' % %& ' ijk dZidtjdtk' m dmM^i 

9 



for i, j, k, I = 1, 2 and 

9ij "j «jj yijjfe "j "j'Jfej &ijkl «j "jkh ^ijklm «j ^jklmi Jijklm ^M^feim- 

Since F^. can be written as the linear combinatorics of tj, ty , t^fe and ty/w, all 
the inner products in ffTTj) and (IT2l can be expressed as linear combinations 

OI yijj 9ijki Qijkh ^ijklj ^ijklm anfl Jijklm- 

Substituting ((TTD and (|I2j) into (|HJ), © and (HDD, and then substituting 
(JHD , © and ( TTU1) into the expression |A(p, r) — | ^- cot(7i(r))df(r), and 
using Maple to conduct all the symbolic calculation, we have 



6 = a /3 = \Ja 2 t 2 (g ll g 2 2 - </? 2 )> & i = °- 
The theorem is proved. □ 

Remark 4.2. The calculation of &o> &i involves a huge number of terms. It 
is almost impossible to finish the derivation by hand. Maple completes all 
the computation in 26 seconds on a PC equipped with a 3.0GHZ Intel(R) 
CPU. The Maple code that conducts all derivation of the theorem is available 
in http : //lsec . cc . ac . cn/~xuzq/maple . html The interested readers are 



encouraged to perform the computation. 

Remark 4.3. It should be pointed out that there is another discrete scheme 

r (4) = 2?r ~ Ei Ji 
' A M (p) ' 

where A M (p) is the area of Voronoi region. Since £\ area(pjpp i+1 ) could be 
approximated by 3 Am (p) under some conditions, for example the conditions 
of Theorem I4.1[ G^ is easily derived from G^ . 

5 A New Discrete Scheme of the Gaussian 
Curvature and Its Convergence 

In this section, we introduce a new discrete scheme for Gaussian curvature 
which converges over the umbilical points and regular vertices with valence 
greater than 4. This is the main result of the paper. We firstly discuss some 
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properties about the discrete mean curvature. Setting ctj 
Si = Zpip i+ ip, we let 



#(i) 



J2i(cotai + cot6i)pPi 



^(cotai + cot 5i)rf 



Zpip^xp and 



(13) 



which is a popular discrete scheme for the mean curvature at vertex p (c.f. 
[T3]). Moreover, the real mean curvature and the real Gaussian curvature at 
p are denoted as H and G respectively. Then, we have 

Lemma 5.1. At the regular vertex p, or the umbilical points, the discrete 
scheme H^ x > converges linearly to the mean curvature H as rj = r\i — > 0. 

Proof. Firstly, let us consider the convergence property at the regular vertex. 
Since p is a regular vertex, cota'+cot^ = 1 + 0(rj 2 ), for any different % and j. 
It follows from equation (jSJ), we have 



Nfcot «j + cot 5i)ppi 



£ 



(cot ati + cot Si)r]fki 



+ 0( V 3 ). 



Hence, 



Secondly, we study the convergence properties at the umbilical points. 
Over the umbilical points, ki = kj = H for any i and j. Hence, 



H m 



Y,i(cot a { + cot Si)ppi 



J2i(cot ai + cot 5i)rjf 
J2i(cot on + cot 5i)rjfki + 0{rf) 



H + 0( V ). 



J^Acotai + cot 5i)r]f 

Combining the two results above, the theorem holds. □ 

Now, we turn to a new discrete scheme for Gaussian curvature. Let 
<Pi '■= Ej=i li and 



G® :-- 



Li\ Ox) 



n 



where 



* == £ 



1 ,Ww-i 



4 sin 7j 2 

COS7i 2 . 2 



1 — COS 2<fii COS 2<y3 i+ i) 



V=£ 



4 
1 



(^ sin2 <ft + ^+i sin2 Vi+i)). 



t-' 4 sin(7i 



Wi+i - 



=^w+*.: 



Then, we have 

Theorem 5.2. For i/ie regular vertices with valence not less than 5, or the 
umbilical points, G^ converges towards the Gaussian curvature G as rji — > 0. 

Proof. We firstly consider the regular vertex case. We set 9{n) := — . Since 
p is a regular vertex, 7j = 9{n) + 0(r] 2 ) for any % according to (jl|). After a 
brief calculation, we have A = A' + 0(r) ), S p = S' p + 0(7/ 4 ), where 



A' 



SL 



16sin(0(n)) 



[2n — ncos(29(n)) — ncos(9(n))]r]' 2 



n 



4sin(0(n)) 



[1 - cos{9(n))]r] 2 . 



Hence, we have 



(2vr - £ 7 , ~ 2 (Sp - A){H^f)/(2A - S p ) 

i 

= (2tt - J> - 2(S' p - A'){H^f)/{2A' - S' p ) + 0(r, 2 ). 

i 

Note that 2ni2i = 1 + 0(rj). According to Theorem 3 in IT], we have 
2vr - J> = A'G + B\k 2 M + k 2 J + o( V 2 ), 



where, B' 



16sin(0(n)) L" ' 2 

Note that S' = A' + 2B' and 



n + § cos(2#(») - f cos(fl(n))]r/ 2 . 



A'G + B'(A& + fc 



U2 ^ 



A'G + B'[(k M + k m ) 2 -2k M k r 
A'G + 45'i/ 2 - 2B'G 
(A' -2B')G + AB'H 2 . 



12 



Hence, 2vr-^. 7i = (A' -2B')G+AB'H 2 + o{r] 2 ). Note that A' = 0{r] 2 ),B' = 
0(r] 2 ) and A' — 2B' ^ provided n ^ 3. According to Lemma 15.11 H^- 1 ' 
converges to the real mean curvature. Hence, we have, when n > 5, 

= 2,-E.7.-4^ 

A'-2B' K ' 

= 2A ,_s> P L + o(l) = G^ + o(l). 

Therefore, G^ converges to the Gaussian curvature. 
Now, let us consider the umbilical point case. For umbilical points, each 
directional is the principal direction. According to Lemma 4 in jl] , we have 

2vr - J2 7i = (AG + (S p - A)k 2 m ) + o( V 2 ) 

i 

over the umbilical points. Since k^ = H 2 = G, we have 
2vt-^ 7j = (AG + (S p - A)k 2 J + o{rf) 

i 

= (AG + 2(S p -A)H 2 -(S p -A)G) + o(ri 2 ). 
Hence, 

G = ^-^, ; 2(S-A)(H^ + q(iH G(5) + o(i) 

The theorem holds. □ 

Remark 5.3. Theorem 15.21 shows that the new scheme G^ 5 ' converges over 
the regular vertex with valence greater than 4. As shown before, the previous 
schemes only converge over the regular vertex with valence 6, and hence the 
new scheme has better convergence properties over the available scheme. 

Remark 5.4. In [8], the authors also prove that the discrete scheme H^ 
converges to the real mean curvature at the regular vertex. However, the 
definition of the regular vertex in [8] is different with our definition. 

Remark 5.5. According to the conclusions above, the Gaussian curvature 
and mean curvature can be approximated over the regular vertex with valence 
greater than 4. Hence, using the formulation k m = H — \J H 2 — G, ku = 
H + ^/H 2 — G, one can approximate the principal curvatures over the regular 
vertex with valence greater than 4. 

13 



6 A Counterexample for the Regular Vertex 
with Valence 4 



In [21], we have constructed a triangular mesh and shown that it is impossible 
to construct a discrete Gaussian curvature scheme which converges for that 
mesh. But the vertex in the mesh is not regular. In this section, we shall 
show that it is also impossible to build a discrete Gaussian curvature scheme 
which converges over the regular vertex with valence 4. 

Suppose the xy plane is triangulated around (0, 0) by choosing 4 points 
qi = (r!,0),q 2 = (0,r 1 ),q 3 = (—r 1 ,0) and q 4 = (0, — ri). For a bivariate 
function f(x, y), the graph of f(x, y), i.e. F(x, y) = [x, y, f(x, y)] T , can be re- 
garded as a parametric surface. Let p = F(0, 0) and p^ = F(qj), i = 1, 2, 3, 4. 
The set of triangles PiPoPi+i forms a triangular mesh approximation of F at 
Po- The triangular mesh is denoted as Mf. When f(x,y) is in the form of 
x 2 + cxy + y 2 where c £ M, it is easy to prove that p := (0, 0, 0) T is a regular 
vertex with valence 4. Moreover, p 1 = (r 1; 0,r 2 ) T , p 2 = (0,r 1 ,r 2 ) T , p 3 = 
(— ri, 0,r 2 ) T , p 4 = (0, — ri,rf) T . Now we show that it is impossible to con- 
struct a discrete scheme for Gaussian curvature which converges over the 
vertex p (See Fig. 2). 

We assume that the discrete scheme for Gaussian curvature involving one- 
ring neighbor vertices of p , which is denoted as G(Mf, p ; Pi, p 2 , P3, P4), is 
convergent for the regular vertex with valence 4 over triangular mesh surface 
Mf, where f(x, y) is in the form of x 2 + cxy + y 2 . 

It is easy to calculate that the Gaussian curvature of F(x, y, z) at p 
is 4 — c 2 . By the convergence property of G(Mf, p ; pi, p 2 , P3, P4) we have 
lim ri ^ G(Mf, p ; p lf p 2 , p 3 , p 4 ) = 4— c 2 . Note that the triangular mesh Mf is 
independent of c, i.e. for any function f(x, y) which is in the form of x 2 +cxy+ 
y 2 , the triangular mesh Mf is the same. Hence, lim ri ^ G(Mf, p Q ; pi,p2, P3, Pi) 
is independent of c. A contradiction occurs. 

Hence, the assumption of G(Mf,p ;pi,p2,P3,Pi) being convergent for the 
triangular mesh Mf does not hold. 



14 



Fig. 2. A sequence of regular vertex with valence n = 4 for the function 
f(x, y) = x 2 + xy + y 2 . At the regular vertex, it is impossible to construct a 
discrete Gaussian curvature scheme which converges to the correct value. 

Remark 6.1. The counterexample in this section justifies the conclusion 
in [TJ , which says that 4 is the only value of valence such that 2ir — £\ ji 
depends upon the principal directions. 

Remark 6.2. An open problem is to find a discrete scheme for Gaussian 
curvature which converges at the regular vertex with valence 3. 



7 Numerical Experiments 

The aim of this section is to exhibit the numerical behaviors of the discrete 
schemes mentioned above. For a real vector a = (a 2 o, an, 002)? we define a 
bivariate function f a (x, y) := a 20 x 2 + a\\xy + a 02 y 2 , and regard the graph of 
the function f a (x, y) as a parametric surface 

F a (x,y) = [x,yJ a (x,y)f ER 3 . 

The Gaussian curvature of F a (x, y) at the origin is 4a 2 octo2 — a n- The domain 
around (0, 0) is triangulated locally by choosing n points: 



q k = h(cos9 k ,sm9 k ), 9 k = 2(k - l)n/n, k = 1, 



n. 



Let pfc = F a (q fc ) and p = (0, 0, 0) T . Hence, the set of triangles {p^PoPfc+i} 
forms a piecewise linear approximation of F a around po- We set e k := 
f a (cos9 k , sin9 k ) and select 



\ 



24 



k> 2 



(14) 



15 



Table 1: The asymptotic maximal error eW(n). 

t-\ \ / \ I 7"oT7 \ I ! a \ / \ I 



n 


eW(n) 


eW{n) 


eW(n) 


e^(n) 


4 


4.6016e + 01 


3.3571e + 01 


3.3570e + 01 


3.3593e + 01 


5 


8.2000e + 00 


9.3792e + 00 


9.3792e + 00 


4.1631e + 01r/ 


6 


1.2226e + 01r/ 


1.2903e + 01r/ 


1.2903e + OI77 


1.1488e + 01r/ 


7 


3.8464e + 00 


4.5783e + 00 


4.5783e + 00 


9.0676e - OI77 


8 


5.8387e + 00 


7.7628e + 00 


7.7628e + 00 


6.5630e + 01r/ 



so that po is a regular vertex. 

We let (jW(po : F a ) denote the approximated Gaussian curvatures of F a 
at po, which is obtained by using the discrete scheme G^\ Suppose A is a 
set consisting of M randomly chosen vectors a. Then, we let 



e®(n) 






|G«(Po:F a )-(4a 20 a 02 



l ii 



)\/M. 



In fact, £^'(n) measures the error of the discrete scheme G^> at the regular 
vertex with valence n. The convergence property and the convergence rate 
are checked by taking l-y = 1/8, 1/16, 1/32, ■ • • . (when k > 2, Z*. can be 
obtained by (fT4")) .) Since p is regular, each edge has the same length 77. Table 
1 shows the asymptotic maximal error e^'(n) for M = 10 4 . Here, the vertex 
valences n are taken to be 4, 5, ■ • ■ ,8. 

From table 1, we can see that all methods work well on valence 6 but 
only new method works well for valence > 5. 

We compute the Gaussian curvature over a randomly triangulated unit 
sphere by the discrete schemes to test their convergent property at the umbili- 
cal points. The vertexes of the random triangulation are uniform distribution 
on the sphere. Fig. 3 shows the random triangulation for the unit sphere. 
Denote the vertices in the random triangulation as p^, i — 1, • • • , N where N 
is the number of the vertices in the random triangulation. We let G^'(pi) 
denote the approximate Gaussian curvature at the vertex p^ which is calcu- 
lated by G( j l Similarly to the above, we use e® = ^ii \G {j) (Pi) - 1)\/N 
to measure the error of discrete scheme G^ and use 77 to denote the average 
length of the edges. Table 2 lists e^' for different N. Moreover, we also use 
e^ to denote the error of the discrete scheme H^ l > for the mean curvature. 
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Table 2: The asymptotic error eW over a sphere with very irregular connec- 
tivity. 



N 


V 


eW 


£ (2) 


£ (4) 


£ ( 5 ) 


e U0 


30 


0.710 


3.798e-01 


1.905e-01 


1.905e-01 


2.126e-01 


2.840e-02 


100 


0.383 


3.517e-01 


5.480e - 02 


5.480e - 02 


1.192e-01 


1.301e-02 


400 


0.196 


2.673e-01 


1.280e - 02 


1.280e - 02 


1.730e-02 


2.600e-03 


1300 


0.109 


2.812e-01 


3.801e-03 


3.801e-03 


6.500e-03 


7.540e-04 


5000 


0.056 


2.669e - 01 


9.648e - 04 


9.648e - 04 


2.703e - 03 


1.893e-04 







v ,' 



Fig. 3. Our test random triangulations. From left to right, the number of 
vertices is 30, 100, 400, 1300, 5000 respectively. 

From these numerical results, we can draw the following conclusions: 

1. For the regular vertices with the valence greater than 4, or the umbilical 
points, the discrete scheme G^ converges to the real Gaussian curvature. 
This agrees with the theoretical result. 

2. At the regular vertices and the umbilical points, the difference between 
G^ 2 ' and G^ A ' is very small. 
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